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Vertically stacked and coupled InAs/GaAs self-assembled quantum dots (SADs) are predicted to 
exhibit a strong non-parabolic dependence of the interband transition energy on the electric field, 
which is not encountered in single SAD structures nor in other types of quantum structures. Our 
study based on an eight-band strain-dependent k ■ p Hamiltonian indicates that this anomalous 
quantum confined Stark effect is caused by the three-dimensional strain field distribution which 
influences drastically the hole states in the stacked SAD structures. 
PACS: 78.67.Hc, 73. 21. La, 71.70.Fk 
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Zero-dimensional semiconductor structure—like InAs 
/GaAs self-assembled quantum dots (SADs)EJ have at- 
tracted— considerable attention because of the new 
physicsoQ of a few electrMi systems and potential appli- 
cations in optoelectronicsO. -Recent experiment on Stark 
effect spectroscopy in SADa£l has demonstrated the exis- 
tence of an inverted electron-hole alignment due to the 
presence of gallium diffusion in InAs SADs, and estab- 
lished a relation between the Stark shift and the vertical 
electron-hole separation. 

The theoretical interpretation of these experimental re- 
sults is based on the assumption that the applied elec- 
tric field can be treated by the second-order perturbation 
theory, which results in a quadratical dependence of the 
transition energy on the applied electric fieldQ, 



E{F) ^E{0)+pF + f3F'^ 



(1) 



where p is the built-in dipole moment and (3 measures 
the polarization of the electron and hole states, i.e., the 
quantum confined Stark effect. While this relation is 
well satisfied in many quantum systems-.ipcluding sin- 
gle SADsQ, and quantum well structuresElB, we show in 
this workptiat it is not valid for vertically coupled SAD 
structurealj where the quantum confined Stark effect de- 
viates significantly from its quadratic dependence on the 
electric field. The reason for this anomalous quantum 
confined Stark effect is due to the three-dimensional (3D) 
strain field distribution in the dots and in the coupling 
region, which controls the localization of hole states in 
the respective SADs, and their sensitivity to external 
field. The existence of this effect is important for basic 
condensed matter physics because it can not be inferred 
from a simple superposition of the electronic properties 
of single SADs. ft is also promising for applications in 
optoelectronics because interband transition energies can 
be significantly modulated by electric fields in quantum 
dot lasers and other photonic devices. 

The insets of Fig. 1 show schematically a single SAD 
structure and a system of two vertically coupled SADs 
that are truncated pyramids separated by a GaAs bar- 



directed from the top to the bottom of the structures. 
Fig. 1 shows the calculated ground state transition ener- 
gies for the single dot and for the stacked structure, as 
functions of electric fields. The electron and hole states 
of the system are obtained from the Schrodinger equaJion 
in the framework of the envelop function formalismEj, 



(Hfc.p -I- \e\Fz)^ = E(t). 



(2) 



Here (f> = {<pi, (/)2, . . . , ^s) is the envelop eigenvector and 
Hfe.p is the k ■ p eight-band Hamiltonia«|-which includes 
the effect of strain and piezoelectricityt3'0. The Hamil- 
tonian is discretized on a three-dimensional grid as a large 
sparse matrix which is solved by Lanczos algorithm. This 
approach has been shown to be reliable, especially in the 
investigation of-the inverted electron-hole alignment in 
SAD structurc£3. 
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FIG. 1. Ground state transition energies as functions of 
vertical electric fields for single dot (top) and stacked double 
dots (bottom). Insets show the dot systems schematically. 
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FIG. 2. Energy levels of stacked double dots, as a function 
of the electric field. The dash-dotted lines are for the single 
dot of the same size. 

The single dot exhibits a nearly perfect quadrati- 
cal dependence on the electric field which is referred 
to the conventional quantum confined Stark eff'ect and 
has been |-p|bserved in many other types of quantum 
structureaJH. The maximum transition energy occurs at 
F « — 50 kV/cm because this homogeneous InAs SAD 
has a positive electron-hole alignment, i.e., the hole is 
at the bottom of the dot for zero electric field. The 
stacked structure shows piecewisely quasi-linear depen- 
dence, with a turning point at a smaller negative elec- 
tric field F « — 22 kV/cm. The lower panel of Fig. 1 
shows that this discrepancy from the conventional quan- 
tum confined Stark effect occurs over the whole range of 
electric field, i.e., for 200 kV/cm < F < 200 kV/cm. In 
addition, the stacked structure exhibits a Stark shift of 
one order of magnitude stronger than the single dot. We 
also notice that the maximum transition energy for both 
systems are roughly identical, because in the stacked 
structure, despite the downshift of the electron states 
due to tunnel coupling, the hole states are characterized 
by an up-shift in energy due to the strain distribution 
(see Fig. 2). Similar simulations have been performed 
for stacked structures with different sizes, different thick- 
nesses of the coupling region, and with or without wet- 
ting layers, where the same kind of aoomalous quantum 
confined Stark effect has been foundE3. 

Fig. 2 shows the energy levels of the stacked struc- 



ture, for both electrons and holes, as a function of electric 
fields. All the energies are given in reference to the top of 
the valence bands of GaAs. For the sake of showing the 
fine structure of the hole spectra near zero electric field, 
we use an energy scale twice as small as for electrons. 
For comparison, the energy levels of the single dot, are 
shown in dash-dotted lines. It is clearly seen that, due to 
the coupling between the stacked structures, the double 
dot system exhibits much richer structure in its energy 
spectrum. 

In the conduction band, the energy spectrum of the 
single dot shows the ground Is state, the two nearly 
degenerate excited 2p states, and the 3d-\ike state with 
weak sensitivity to the electric field. The coupled dot sys- 
tem results in bonding and antibonding states originating 
from these states, that in the absence of electric field, are 
identified as ls+, ls~, 2p+, 2p~, 2s~^,3d~^, ■ ■ ■ where +(— ) 
denotes bonding (antibonding) states. Except for the 
ground state, all the excited states are seen to have cross- 
ings or anticrossings with other states, which reorders the 
states at high fields. In the valence bands of the single 
dot, the four hole levels {hf, /i|, /i| and /i|) are seen to 
show quadratical dependence on the electric field. 

The most dramatic feature in the valence bands of the 
stacked structure is the quasi-linear dependence of the 
hole levels on the electric field for intermediate and strong 
field intensities, especially the ground hole state hf. In 
addition, we could not find any bonding or antibonding 
hole states for these structure as the ground hole state hf 
and first excited state /i2 are localized in the bottom dot 
and the top one, respectively. However, at small electric 
fields, the hole levels show small fluctuations due to mu- 
tual anticrossings. It is also noticed that the magnitude 
of the Stark shift for holes in the stacked structure is sig- 
nificantly larger than in the single dot, which is respon- 
sible for the large shift in the transition energy shown in 
Fig. 1. 

The main reason for the hole states to behave so dif- 
ferent in the stacked structure than in the single dot is 
traced in the 3D strain field distribution. Fig. 3(a) shows 
the profiles and the corresponding contour plots (insets) 
for the hydrostatic (H) and the biaxial (B) components 
of the strain field, which are defined as 



H 



B'- - (e. 



,? + {^y 



yj -\- y^zz ^xx) 



(3) 



Fig. 3(b) shows the diagrams of thCj-cpuduction and va- 
lence band edges which are given byll3 



AUc = AEc + ac- H 
AUy = AEy -a^- H -b- B/V^ 



(4) 



where AEc (AEy) is the conduction (valence) band off- 
set, Oc (at,) is the conduction (valence) band hydrostatic 
deformation potential parameter, b is the valence band 
shear deformation potential parameter. Hence, while 
electrons are only sensitive to hydrostatic strain, holes 
are also and mostly affected by biaxial strain. 
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FIG. 3. (a) Biaxial (upper curve) and hydrostatic (lower 
curve) components of the strain field in the single (dash lines) 
and stacked (solid lines) SAD structures through the center 
of the dots along the growth direction (z). The zero on the 
horizontal axis is fixed at the top of the upper dot. (b) Band 
diagrams for the single (dash lines) and stacked (solid lines) 
SAD structures. Inset: Schematic view of the ground hole 
state confined by a triangular potential in the bottom dot 
at zero electric field (solid lines) and a positive electric field 
(dash lines). 

In the stacked structure, the hydrostatic negative 
strain field resides entirely inside the dots, and is a little 
stronger than in the single dot, while it almost vanishes 
in the coupling region (lower curve in Fig. 3a). Therefore, 
the coupling region is seen by electrons as a conventional 
tunneling barrier, which results in bonding and antibond- 
ing states as shown in Fig. 2. The biaxial positive strain 
field is however seen very differently in the stacked struc- 
ture from that in the single dot (upper curve in Fig. 3a). 
First, it is smaller than in the single dot, leading to hole 
levels with higher energies in the stacked system (see 
Fig. 1). Second, unlike the hydrostatic strain, the biaxial 
strain retains a substantial value in the coupling region 
in the stacked system (top inset of Fig. 3a), which no- 
ticeably reduces the barrier height in the valence bands. 
Third, the biaxial strain profile in the stacked structure 
is inverted in the upper dot compared to the lower one 
although symmetric with respect to a median plane be- 
tween the dots. 




FIG. 4. Probability density isosurfaces of ground states of 
electron (top panel) and hole (two bottom panels) at differ- 
ent electric fields. From left to right, the electric field is (a) 
F = -200, 0, and -)-200 kV/cm; (b) F = -21.8, -21.6, and 
— 21.4 kV/cm, respectively. 

The most important point is that the biaxial strain 
defines two triangular confining potentials in each dot of 
the stacked structure, while in the single dot the valence 
band edge profile exhibits a much smoother slope (see 
Fig. 3b). It is seen that the strain- induced triangular 
potential is inverted with respect to the base in the two 
dots: In the bottom dot, the valence band edge is higher 
at the base than at the top, while it is opposite in the top 
dot. Therefore, the favorable combination of the lower 
strain-induced potential and the wider base result in a 
weaker hole confinement in the bottom dot than in the 
top dot, thereby localizing the hole ground state in the 
lower dot. The inset of Fig. 3(b) schematically shows 
the potential in the lower dot of the stacked structure 
together with the ground hole state at different electric 
fields. At zero electric field, the strain-induced built-in 
potential variation in the dot is 34 meV and the ground 
hole state is localized approximately one fourth of the 
lower dot height from the bottom base. 

This situation is illustrated in the middle panel of 
Fig. 4(a) where we plot the probability density isosur- 
faces of ground electron and hole states at dilferent elec- 
tric fields. In contrast, it is seen that the ground electron 
state is the bonding ls+ state extending quasi-equally 
in both dots and in the coupling region. Under strong 
positive electric field, the ground electron state is seen 
in the right panel of Fig. 4(a) to become a Is-like state 
entirely localized in the upper dot. However, the ground 
hole state is found to behave very differently. Because the 
hole state is localized inside this triangular potential, as 
schematically shown in the inset of Fig. 3(b), its energy 
level changes approximately by the same amount as ex- 
ternal variation in the dot, i.e^he potential drop arising 
from the external electric fieldtll. For example, the varia- 
tion of ground hole state energy between F = 200 kV/cm 
and zero field is about 70 meV, which is about the same 
value as the potential drop over the bottom dot caused 
by electric field F — 200 kV/cm. Therefore it is seen 
in Fig. 2 that the energy of the hole state has an al- 
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most linear dependence on the electric field. This also 
explains the quasi-linear dependence of the ground state 
transition energy as shown in Fig. 1. 

In Fig. 3(b), the ratio between the tunneling barrier 
and band offset in the valence bands is much smaller than 
that in the conduction bands. Consequently, it is much 
easier for holes to tunnel through the barrier at negative 
electric fields than electrons. Fig. 4(b) illustrates the 
transition of the ground hole state from the bottom dot 
to the top one occuring at a small negative electric field 
F ~ —21.6 kV/cm. As the biaxial strain field has dif- 
ferent distribution in the respective dots, the hole states 
exhibit different probability density profiles when local- 
izing in different dots. In Fig. 2, this transition is seen 
as an anticrossing between the ground hole state and the 
first excited state, which also accounts for the anomalous 
quantum confined Stark effect. 

Although the anomalous quantum confined Stark effect 
invalidates Eq. (1) over the whole range of electric field, it 
is possible to write two separate equations that describe 
the dependence of the transition energy on electric field 
for the stacked structure, except for the electric fields at 
which anticrossings occur. 



E-{F) 



EiOj+p+F + P+F"^, 
= E' +p-F + (3-F'\ 



(5) 



where E' is a fitting parameter, E^{F) is for F > 
-21 A kV/cm, and E-{F) is for F < -21.8 kV/cm. 

For the right branch, p^/|e| = —3.65 nm is the built- 
in dipole moment of the stacked structure at zero field. 
From the middle panel of Fig. 4(a), the ground electron 
and hole states are seen to be separated by a distance 
larger than the coupling region (1.8 nm), and much larger 
than the electron-hole separation in the single dot, which 
results in a dipole moment of more than one order of 
magnitude larger than that in the single dot. For the 
left branch, p~ « —p'^ nm has a positive sign because 
of the inverted electron-hole alignmentllJ in the stacked 
structure at negative electric fields. 

The quantum confined Stark effect coefficient is defined 
hy (3 = Pe — Ph where f3e {(ih) can be related to oscillator 
strength of opticjal intraband transitions in conduction 
(valence) bandsQllj. For instance, (3e is given by 



/3e = -;3^E/i-«/(i?n-i?i)', 



2mo 



(6) 



n>l 



where toq is the bare electron mass, fi^n is the oscilla- 
tor strength for the intraband transition from the ground 
state to the n-th state, with polarization along z direc- 
tion. Ph has a similar expression. Single SAD structures 
have been showp_.to have very weak z-polarized intra- 
band transitionsES, therefore both /3e and Pu are small. 
In the stacked structure, the strengths of the z-polarized 
intraband transitions in the valence bands are similar to 
those in the single dot. i.e., Ph is small, which could be 
observed in Fig. 2 where the hole state energies show a 



quasi-linear dependence on the electric field. However, 
in the conduction bands of the stacked structure, there 
are several strong intpiband transitions, especially the 
Is"*" — > ls~ transitionllS. This results in a much larger 
/3e than in single dots, and explains the magnitude of 
the quantum confined Stark effect coefficients in both 
branches (/3^/|ep « — 80 nxs? /eV) which are six times 
larger than in single dots, and are responsible for the 
'bowing' of the Stark shift in Fig. 1. 

In conclusion, we have shown that vertically stacked 
InAs/GaAs self-assembled quantum dots exhibits, in ad- 
dition to a much larger Stark shift than single dots, a 
strong non-parabolic dependence of interband transition 
energy on the electric field, not been encountered in other 
types of quantum structures. We have demonstrated that 
the 3D distribution of the biaxial strain field is mainly 
responsible for this anomalous quantum confined Stark 
effect. 
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